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We present STAR measurements of the azimuthal anisotropy parameter vi for pions, kaons, 
protons, A, A, 3 + S, and Q + Q, along with Vi for pions, kaons, protons, and A + A at mid-rapidity 
for Au+Au collisions at ^/s NN = 62.4 and 200 GeV. The V2(pr) values for all hadron species at 62.4 
GeV are similar to those observed in 130 and 200 GeV collisions. For observed kinematic ranges, V2 
values at 62.4, 130, and 200 GeV are as little as 10%-15% larger than those in Pb+Pb collisions at 
y/s NN = 17.3 GeV. At intermediate transverse momentum (pr from 1.5-5 GeV/c), the 62.4 GeV 
V2(pt) and v^pt) values are consistent with the quark-number scaling first observed at 200 GeV. 
A four-particle cumulant analysis is used to assess the non-flow contributions to pions and protons 
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and some indications are found for a smaller non-flow contribution to protons than pions. Baryon vi 
is larger than anti-baryon V2 at 62.4 and 200 GeV perhaps indicating either that the initial spatial 
net-baryon distribution is anisotropic, that the mechanism leading to transport of baryon number 
from beam- to mid- rapidity enhances V2, or that anti-baryon and baryon annihilation is larger in 
the in-plane direction. 

PACS numbers: 25.75.Ld, 25.75.Dw 



I. INTRODUCTION 

In non-central heavy- ion collisions, the overlapping 
area has a long axis and a short axis. Re-scattering 
amongst the system's constituents converts the ini- 
tial coordinate-space anisotropy to a momentum-space 
anisotropy [3, & HI- The spatial anisotropy decreases 
as the evolution progresses so that the momentum 
anisotropy is most sensitive to the early phase of the evo- 
lution — before the spatial asymmetry is washed-out [3| • 

Ultra-relativistic Au+Au collisions at Brookhaven 
National Laboratory's Relativistic Heavy Ion Collider 
(RHIC) H are studied in part to deduce whether quarks 
and gluons become deconfined during the early, high 
energy-density phase of these collisions. Since the az- 
imuthal momentum-space anisotropy of particle produc- 
tion is sensitive to the early phase of the collision's 
evolution, observables measuring this anisotropy are es- 
pecially interesting. The azimuth angle (</>) depen- 
dence of particle momentum distributions can be ex- 
pressed in the form of a Fourier series: dN/d(j) oc 1 + 

2v n cos n ((/) — ^rp), where ^pp is the reaction- 
plane angle 0,0- The Fourier coefficients v n can be mea- 
sured and used to characterize the azimuthal anisotropy 
of particle production. 

Measurements at two higher RHIC energies (y/s NN — 
130 and 200 GeV) established that charged hadron V2 
rises with px for pp < 2 GeV/c and then saturates la Bl- 
As predicted by the hydrodynamic calculations [HJ, LL1( 
- where local thermal equilibrium is assumed — V2 
at low pt (pt < 1 GeV/c) shows a characteristic de- 
pendence on particle mass [HI, [l3j]. The V2 values at 
y/s NN = 130 and 200 GeV are as large as those predicted 
by hydrodynamic calculations. The V2 values measured 
at y/ s NN = 17.3 GeV [3], the top energy of the Super 
Proton Synchrotron (SPS) at CERN, however, are be- 
low the hydrodynamic models predictions. In this paper, 
we compare V2 at ^s NN — 17.3 and 200 GeV to new 
measurements at y^^T = 62.4 GeV that provide a link 
between the top RHIC energy and the top SPS energy. 

In 200 GeV collisions, kaon, proton, A + A and E + 5 
V2(pt) at intermediate pp depends on the number n q 
of constituent quarks in the corresponding hadron [i~5l |. 
A scaling law — motivated by constituent-quark coales- 
cence or recombination models — can account for the 
observed splitting between bar yon s and mesons for i>2 
in this intermediate region [la, [l6| . Within these mod- 
els, hadron w 2 (^2) is related to the i>2 of quarks (vf) 
in a quark-gluon phase by the relationship: ^(Pt) ~ 



n q v\(n q p T ) [13] • Intermediate pp baryon yields also in- 
crease with collision centrality more rapidly than meson 
yields [HI, (HI : a behavior also expected from coalescence 
or recombination models [l7} ■ These models suggest that 
the large V2 values at intermediate pp are developed dur- 
ing a pre-hadronic phase — a conclusion supported by 
the recent discovery that multi-strange bar yon s, thought 
to have smaller hadronic cross-sections [THHo)], attain i> 2 
values apparently as large as protons or hyperons [2lj . 
Measurements of V2 for identified particles may, there- 
fore, help reveal whether V2 is developed in a deconfined 
quark-gluon phase and can test whether these possible 
conclusions may still be valid at lower ^s NN . 

This article is organized as follows: in Sec. [H] we 
briefly describe the STAR detector. The analysis pro- 
cedures are presented in Sec. IIIII In Sec. IIVI we present 
our results. This section includes subsections discussing 
systematic uncertainties, baryon versus anti-baryon u 2 , 
quark- number scaling, the fourth harmonic i>4, and the 
collision energy dependence of V2- Our conclusions are 
then presented in Sec. [Vj 



II. EXPERIMENT 

Our data were collected from Au + Au collisions at 
= 62.4 and 200 GeV with the STAR detector [13]. 
STAR's main time projection chamber (TPC) [23[ was 
used for particle tracking and identification with supple- 
mentary particle identification provided by time-of-fiight 
detectors (TOF) [24[. We analyzed events from a cen- 
trality interval corresponding to 0%-80% of the hadronic 
interaction cross-section. As in previous STAR analy- 
ses [lH, we define the centrality of an event from the 
number of charged tracks in the TPC having pseudo- 
rapidity \ri\ < 0.5, pt > 0.2 GeV/c, a distance of clos- 
est approach to the primary vertex (DCA) less than 2 
cm, and more than 10 measured space points [25j |. Only 
events with primary vertices within 30 cm of the TPC 
center in the beam direction were analyzed. 

STAR's main TPC covers the approximate pseudo- 
rapidity region |?7| < 1.2 (for collisions at its center) and 
2ir in azimuth angle. A 0.5 Tesla magnetic field allows 
charged particle pt to be measured above 0.1 GeV/c. 
At the time of data taking the TOF detectors covered 
— 1 < i] < and 7r/15 in azimuth angle. Their timing 
resolutions are ~ 110 ps so that pions and kaons can be 
distinguished for pp < 1-8 GeV/c and protons can be 
identified up to pp = 3.0 GeV/c. 
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III. ANALYSIS 

We identify particles using three different methods: 
measurement of specific ionization-energy-loss per unit 
length in the TPC gas (dE/dx), time-of- flight measure- 
ments, and weak-decay vertex finding. dE/dx measure- 
ments for a particle with a given momentum are used for 
identification at low pt and in the relativistic-rise region 
(jPt > 2.0 GeV/c) where dE/dx increases logarithmically 
with 07 (see Ref. ^ and Fig. 26 in Ref. [27]). The pion 
sample in the relativistic-rise region is selected based on 
the deviation between the measured dE/dx of each track 
and the expected dE/dx for a pion in units of Gaussian 
standard deviations (nov). For pt > 2.0 GeV/c, pions 
are selected with nov > (the top half of the distribu- 
tion). In this case the purity is estimated to be 98%. 

The V2 of protons is measured in this region by fitting 
the dE/dx distribution with peaks centered at the pre- 
dicted dE/dx values. From these fits we can derive the 
relative fractions of pions (f w ), kaons (fa) and protons 
(f p ) as a function of dE/dx. We then measure V2 for all 
tracks and plot it versus the dE/dx of the track. Once 
the relative fractions of each particle are known for each 
value of dE/dx, and vi is know as a function of dE/dx, 
V2(dE / dx) can be fit with function: 

V 2 (dE/dx) = f„V 2 ,-n + fKV2,K + fpV2,p, (1) 

where the v-i values for each species (w2,ir, V2,k, and V2, p ) 
are parameters in the fit and fx, and f p which are 
extracted from the dE/dx distribution, are part of the 
fit function. In the relativistic rise region, kaons do not 
dominate the dE/dx distribution for any value of dE/dx, 
so their d 2 values are poorly constrained and are not pre- 
sented here. We estimated the systematic error on the 
proton V2 by varying the relative fractions of the different 
particles within reasonable limits. The relative change in 
the proton V2 (SV2/V2) was less than 3%. The shape and 
width of the peaks are determined from samples of par- 
ticles identified by other means, e.g. TOF measurements 
and Kg or A decay daughters. 

The reaction-plane direction is estimated for each 
event from the azimuthal distribution of charged tracks. 
We select tracks using criteria similar to those in Ref. || . 
To avoid self-correlations, we subtract the contribution 
of a given particle from the total reaction-plane vector. 
For particles identified through their decays, we sub- 
tract the contributions of all the decay products. The 
reaction-plane resolution is estimated using the sub-event 
method [28| and we correct the observed v 2 to account 
for the dilution caused by imperfect resolution. The res- 
olution depends on the number of tracks used in the cal- 
culation and the magnitude of V2 , and therefore depends 
on centrality. The resolution for y/s~^/ — 62.4 GeV col- 
lisions is reduced relative to y/s NN = 200 GeV collisions 
by w 30%. For 62.4 GeV Au+Au collisions it reaches a 
maximum value of approximately 0.73 in the 10%— 40% 
centrality interval. 



IV. RESULTS 

In Fig.rfjthe minimum-bias, mid-rapidity V2 values are 
shown for inclusive charged hadrons, pions, kaons, pro- 
tons, A + A, 3 + S, and f2 + fi. The gross features of 
V2 at y / s n ^7 = 62.4 GeV are similar to those observed 
at y/sf^ = 200 GeV pj, For p T < 1.5 GeV/c, a 
mass hierarchy is observed with V2 smaller for heavier 
particles. The px and mass dependencies are qualita- 
tively (not necessarily quantitatively) consistent with ex- 
pectations from hydrodynamic calculations that assume 
the mean- free- path between interactions is zero (lp| . For 
Pt > 2 GeV/c, V2 reaches a maximum, the mass order- 
ing is broken, and «2 for protons and hyperons tend to be 
larger than for either pions or kaons. The V2 values for 
protons and A + A above pr — 2 GeV/c are similar. In 
this region, the multi-strange baryons also exhibit t>2 val- 
ues similar to protons. While hadrons containing strange 
quarks are expected to be less sensitive to the hadronic 
stage, we do not see a statistically significant reduction 
in the V2 values of strange baryons compared to protons. 
Statistical uncertainties, however, still do not exclude the 
possibility of some strangeness content dependence for 
V2- If V2 or its hadron species dependence is developed 
through hadronic interactions, V2 should depend on the 
cross-sections of the interacting hadrons (with hadrons 
with smaller cross-sections developing less anisotropy). 
The large «2 values for 3 + E and f2 + f2 are consistent 
with V2 having been developed before hadronization. 

The centrality dependence of identified hadron V2{pt) 
for ^fs NN = 62 A GeV is shown in Fig. [2l Similar pr and 
mass dependencies are observed for each of the centrality 
intervals: 0%-10%, 10%-40%, and 40%-80%. The data 
from the 0%-10% interval are most affected by non-flow 
effects Q while the 10%-40% interval is least affected 
by these uncertainties. The particle-type dependence of 
non-flow will be discussed in the following section. 

A. Systematic Uncertainties 

Systematic uncertainties are shown in Fig. [1] as bands 
around V2 — 0. The errors are asymmetric. The portions 
of the band above zero represent the negative errors so 
that the difference between the measurement and zero 
is more visually evident. These uncertainties take into 
account effects from weak-decay feed-down, tracking ar- 
tifacts, detector artifacts, and non-flow effects. Non-flow 
effects are dominant. In Fig. [2 the tracking and non- 
flow systematic uncertainties are shown as bands around 
v 2 = and the weak-decay feed-down uncertainties are 
included in the error bars on the pion data points. 

The number of tracks coming from weak-decays that 
are included in the V2 analysis depends on the experimen- 
tal setup and track selection criteria. Pions produced in 
Kg, A, or A decays tend to be distributed at low pt 
with V2 values larger than the pions from other sources. 
We have calculated their effect on the observed pion V2- 
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FIG. 1: (color online). Minimum-bias (0-80% of the collision cross section) «2(pr) for identified hadrons at |?j| < 1 from 
Au+Au collisions at ^/s NJV = 62.4 GeV. To facilitate comparisons between panels, V2 values for inclusive charged hadrons are 
displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the 
identified particles are shown as shaded bands around vi = 0. 



We assume exponential mr spectra for Kg and A with 
inverse slope parameters of 285 and 300 MeV respec- 
tively. For relative abundances, we take K s /(n + + ir~) 
and (A + ~K)/(tt + + n~) ratios of 0.06 and 0.054 respec- 
tively. The V2 of Kg and A are taken from data. We then 
use a full detector simulation to estimate what fraction 
of the weak-decay products will fall within our detec- 
tor acceptance and pass our track selection criteria. We 
find that for our analysis, feed-down will increase Vi by 
approximately 13% (as a fraction of the original V2) at 
Pt = 0.15 GeV/c. The increase falls to approximately 
3% relative at px — 0.25 GeV/c and is negligible for 
Pt > 0.4 GeV/c. Modifications to the observed proton 
V2 from A and A decays are neglegible due to the simi- 
larity of proton and hyperon v-i- 

t>2 measurements can also be distorted by anti- 
correlations that arise from tracking errors (e.g. track- 
merging and hit-sharing) . These anti-correlations can be 
eliminated by correlating tracks with r\ > (77 < 0) with 



an event plane determined from tracks at r) < —0.15 
(57 > 0.15) (77-subevents). This method also has a differ- 
ent sensitivity to the spurious correlations arising from 
jets and resonance decays (non-flow effects discussed in 
the next paragraph). In this paper, ?7-subevents are used 
to analyze pion, Kg, proton and A + A V2- The remain- 
ing systematic uncertainties from detector artifacts are 
estimated by comparing data taken with different field 
settings: 0.5 Tesla (full-field) and 0.25 Tesla (half-field). 
The STAR experiment did not collect half-field data dur- 
ing the 62.4 GeV data taking period so we use the 200 
GeV data to estimate the uncertainties in the 62.4 GeV 
measurements. From these studies, we assign an uncer- 
tainty to V2 for all particles of ±0.0035 (absolute). 

The dominant systematic uncertainties in V2 measure- 
ments arise from correlations unrelated to the reaction 
plane (thought to be primarily from correlations between 
particles coming from jets or resonance decays or other 
correlations intrinsic to p+p collisions) . When t>2 is mea- 
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FIG. 2: (color online). The unidentified charged hadron, charged pion, K s , charged kaon, proton and A + A V2 as a function of 
p T for 10%-40%, 0%-10% and 40%-80% of the Au+Au interaction cross section at ^/s^7 = 62.4 GeV. Weak-decay feed-down 
errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands 
around V2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above 
(below) zero represents the negative (positive) error. 



sured using an event-plane analysis (v2{EP}) [281, these 
correlations can bias the experimental estimation of the 
reaction plane (the event-plane angle) and change the 
apparent t>2 values (non-flow). A four-particle cumu- 
lant analysis of vi (u2{4}) is less sensitive to non-flow 
effects than a standard analysis but yields larger statisti- 
cal uncertainties 0. Although «2{4} has been shown to 
significantly reduce non-flow uncertainties, some sources 
of uncertainty may remain: e.g. if V2 fluctuates from 
event-to-event i>2{4} may yield values smaller than the 
mean vi [29| . The magnitude of possible non-flow cor- 
relations for unidentified charged hadrons is discussed in 
Refs. [1, [l3j]. Here we also discuss variations of non-flow 
effects between different hadron types. 

For pt < 1 GeV/c, a four-particle cumulant analysis is 
carried out for pions and protons identified with greater 
than 98% purity. To study the hadron-type dependence 
of non-flow effects at intermediate px, we analyze two 
samples of charged hadrons at 2.4 < pt < 3.6 GeV/c: 
one with nov > 0, the other with —5 < nov < —2.5. 



Data from the 10%-40% centrality interval are used. For 
na„ > 0, approximately 98% of the charged tracks are 
pions. For —5 < nov < —2.5, the sample contained ap- 
proximately 75% protons, 19% kaons and 6% pions. The 
ratio of the event-plane t>2 (^{EP}) to the cumulant vi 
(t>2{4}) for the pion sample and the proton sample are 
listed in Table [0 In the Pt region below 1 GeV/c, pro- 
ton V2 does not appear to manifest any non-flow correla- 
tions for either energy. For pions in this region, however, 
non-flow correlations seem to account for 10% of the V2 
measured with the event-plane analysis. 

At intermediate px, V2{EP}/v2{4:} is greater than 
unity for protons and pions. This shows that non-flow 
correlations increase the observed V2{EP} for both pro- 
tons and pions. At 62.4 GeV, the increase is the same 
(within errors) for both particles. With the larger 200 
GeV data set however, we observe a larger non-flow frac- 
tion for pions than protons: the pion V2{E P} / U2{4} = 
1.22 ± 0.02 and v 2 {EP} /v 2 {^} for the proton sample 
= 1.16 ± 0.02. Pion V2{EP}, therefore, appears to be 
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more susceptible to non-flow correlations than vi for par- 
ticles in the proton sample. 



TABLE I: The ratio V2{EP}/« 2 {4} (v 2 from a standard event- 
plane analysis over V2 from a four-particle cumulant analysis) 
for the centrality interval 10%-40% in three pr ranges (units 
for pt are GeV/c). The sample from 2.4 < pr < 3.6 GeV 
labeled as protons contains contamination from pions (6%) 
and kaons (19%). 
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B. Baryon vs. Anti-baryon v 2 

To our knowledge, no prediction for a difference be- 
tween baryon and anti-baryon i> 2 exists in the literature. 
Previous measurements at RHIC of identified baryon 
V2 reported no differences between A and A vi or be- 
tween proton and anti- proton v 2 . Typically the particle 
and anti-particle samples were combined. These mea- 
surements were made with smaller data samples and at 
higher energies where the anti-baryon to baryon yield 
ratios are much closer to unity. Several scenarios can 
lead to a difference between anti-baryon and baryon i> 2 
that is larger when the anti-baryon to baryon yield ratio 
is smaller: (1) baryons may develop larger momentum- 
space anisotropies through multiple rescattering as they 
are transported to mid-rapidity, (2) if the initial spatial 
net-baryon density is anisotropic, flow developing in a 
later stage could convert that spatial anisotropy to an 
observable momentum-space anisotropy, and (3) anni- 
hilation of anti-baryons in the medium can reduce the 
anti-baryon yield, with the reduction larger in the more 
dense, in-plane direction than the out-of-plane direction. 
We consider scenario (1) and (2) to be distinct. In sce- 
nario (1), extra V2 is built up while the baryons are being 
transported to mid-rapidity while in scenario (2) the i> 2 
is established through rescattering after the baryons are 
transported to mid-rapidity. 

In Fig. [3] we show the ratio of A v 2 to A u 2 ■ The data 
are from minimum bias Au±Au collisions at 62.4 and 
200 GeV. The bands on the figure represent the mean 
values of the ratios which are respectively 0.948 ± 0.014 
and 0.971 ±0.005 for 62.4 and 200 GeV. In the measured 
range, the A v% is systematically smaller than the A V2 for 
both energies and within errors is approximately px inde- 
pendent (fitting the data with a straight line yields slopes 
of 0.014± 0.028 (GeV/c)- 1 and 0.017± 0.010 (GeV/c)- 1 
respectively for 62.4 and 200 GeV data). The difference 
between A and A V2 is larger at 62.4 GeV, where the 
A to A yield ratio is smaller. Taking into account the 
A/A yield ratios (measured to be 0.532 ± 0.014 at 62.4 
GeV and 0.77 ± 0.05 at 200 GeV/c (H), we find that at 



FIG. 3: (color online). The ratio of A v 2 to A v 2 . The 
data are from minimum bias Au±Au collisions at -js NN = 
62.4 and 200 GeV. The bands show the average values of the 
ratios within the indicated pr ranges. 



62.4 GeV the net A V2 (the asymmetry of the quantity 
A -_A) is 12% ± 3% larger than the v 2 of all other As 
or As. At 200 GeV it is 13% ± 4% larger. The larger A 
i>2 is not anticipated from the RQMD hadronic transport 
model [3l| where at mid-rapidity, the ratio of anti-proton 
v 2 to proton t> 2 is 1.148 ± 0.084 and the ratio of A V2 to 
A V2 is 1.142 ± 0.123. We note however that this model 
does not reproduce the overall magnitude of V2 at this 
energy either. 
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FIG. 4: (color online). The pr integrated ratio of A v 2 to 
A v 2 for three centrality intervals: 0%-10%, 10%-40%, and 
40%-80%. The data are from Au±Au collisions at ^s NN = 
62.4 and 200 GeV. 

In Fig. H] we display the centrality dependence of the 
Pt integrated A w 2 to A w 2 ratio. A Monte-Carlo Glauber 
model is used to convert the centrality intervals de- 
fined by multiplicity into mean impact parameter val- 
ues. Given the errors we are unable to make a definitive 
statement about a possible dependence of the ratio on 
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centrality. 



C. Quark- number Scaling 

Models of hadron formation by coalescence or recombi- 
nation of quarks successfully reproduce many features of 
hadron production in the intermediate pt region (1.5 < 

p T < 5 GeV/c) HE ES G3- These models nnd that 
at intermediate px, «2 may follow a quark- number (n q ) 
scaling with V2{pt / n q ) / n q for most hadrons falling ap- 
proximately on one curve. In these models, this uni- 
versal curve represents the momentum-space anisotropy 
developed by quarks prior to hadron formation. This 
scaling behavior was observed in Au+Au collisions at 
200 GeV 15j. Approximate quark number scaling of v 2 



also exhists in RQMD models where the scaling is re- 
lated to the additive quark hypothesis for hadronic cross- 
sections Hl|- The RQMD model, however, under- 
predicts the value of i>2 by approximately a factor of two. 
Pre-hadronic interactions are therefore thought necessary 
to generate a V2 as large as that observed at RHIC. If V2 
is predominantly established in this pre-hadronic phase, 
the hadronic cross-sections might not play a dominant 
role in establishing the particle- type dependence of v% . 

Fig.[5]shows V2 scaled by the number of valence quarks 
in the hadron (n q ) as a function of pT/n q (a) and (tot — 
m )/n q (b) for identified hadrons at y / s^" = 62.4 GeV. 
A polynomial function has been fit to the scaled values 
of v 2 for all particles except pions, which, for reasons dis- 
cussed below, may violate the scaling. To investigate the 
quality of agreement between hadron species, the data 
from the top panel are scaled by the fitted polynomial 
function and plotted in the bottom panels (c) and (d) of 
Fig. In panel (c), for pT/n q > 0.6 GeV/c, the scaled 
v% of Kg, K ± , p+p and A + A, lie on a single curve, 
within errors. The 62.4 GeV data for these species are 
therefore consistent with the scaling observed in 200 GeV 
collisions. The 3 + 3 V2 may lie below the curve but the 
current errors do not permit a strong conclusion regard- 
ing deviations between 3 + 3 and p + p or A + A V2 . At 
Pr/n-q < 0.6 GeV/c, the scaling breaks down. 

It was shown that for 200 GeV at — mo < 
0.8 GeV/c 2 , V2(ttit — mo) is a linear function and in- 
dependent of hadron mass [33[ . In Fig. [5] panels (b) and 
(d) we combine tot scaling and n q scaling so that a single 
curve can be used to approximately describe V2 through- 
out the measured range. This is the same scaling as used 
in Ref. (34[ where the figures are labelled KEt (tot — mo 
is the transverse kinetic energy). This combined scaling 
works because in the range where «2 is a linear function 
of tot — mo, dividing by n q does not alter the shape of 
the curve. Once it is observed that V2 for all particles 
follow the same linear function for mx — mo, the scaling 
of V2(mT — mo) with n q becomes trivial. At higher pt, 
i>2 is only weakly dependent on px so that changing the 



axis variable from pt / n q to (tot - 
the scaling significantly. 



mo)/n q does not effect 



Pion V2 deviates significantly from the fit function in 
both panels (a) and (b). The contribution of pions from 
resonance decays to the observed pion v 2 may account 
for much of the deviation for px < 1.5 GeV/c [Hj]. For 
Pt > 1.5 GeV/c, non-flow correlations discussed previ- 
ously may contribute to the deviation. From the results 
in Table U we conclude that non-flow effects tend to be 
larger for pions than protons. Particularly for the 200 
GeV data, removing non-flow contributions will increase 
the difference between pion and proton V2 and will im- 
prove the agreement between pion V2/n q and V2/n q for 
the other measured particles. It has also been suggested 
that constitucnt-quark-number scaling may be violated 
for pions because the pion mass is much smaller than the 
masses of its constituent-quarks. This implies a larger 
binding energy and a wider wave-function for the pion. 
As a result, the approximation that hadrons coalesce 
from constituent-quarks with nearly co-linear momenta 
is broken [35| . 

Fig. [6] shows V2/n q versus (tot — mo)/n q for 0%-10%, 
10%-40%, and 40%-80% most central Au+Au collisions 
at y/SNN = 62.4 GeV. V2/n q for each centrality interval 
is scaled by the mean eccentricity of the initial overlap re- 
gion. The eccentricity is calculated from the mean x and 
y positions of the participating nucleons using a Monte- 
Carlo Glauber model. The coordinate system is shifted 
and rotated so that (0, 0) is located at the center-of-mass 
of the participants and the eccentricity is the maximum 
possible. This is referred to as the participant eccentric- 
ity (£p a rt)- Since the true reaction plane is not know, our 
V2 measurements are sensitive to e par t [3(| ■ For the 0%- 
10%, 10%-40%, and 40%-80% centrality intervals the 
(e P ar t) values respectively are 0.080, 0.247, and 0.547. 

The tot — too and n q scalings shown for minimum bias 
data in Fig. [5] are also valid within the specific centrality 
intervals shown in Fig. [6l Early hydrodynamic calcula- 
tions predicted that V2 should approximately scale with 
the initial spatial eccentricity of the collision system [ij . 
V2/{£ P art) contradicts these expectations and rises mono- 
tonically as the centrality changes from peripheral to 
central. This indicates that central collisions are more 
efficient at converting spatial anisotropy to momentum- 
space anisotropy. 



D. Fourth Harmonic V4 

Higher order anisotropy parameters (t>4, vq, etc.) may 
be sensitive probes of hydrodynamic behavior and the 
initial conditions of the collision system [37| ■ The authors 
of Ref. [38| argue that values of the ratio v^jv\ larger than 
0.5 indicate deviations from ideal fluid behavior. When 
measured for identified particles, higher harmonics can 
also test quark-number scaling [39( . v 4 and vq for charged 
hadrons at 200 GeV are shown in Ref. 14011 . Identified 
particle D4 at 200 GeV is shown in Ref. [13 |. In Fig. [7] 
(top panels) we plot pion, kaon, anti-proton and A + 
A V4 for y 7 ^^" = 62.4 GeV, where the standard event- 
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FIG. 5: (color online). Identified particle i>2 from minimum bias collisions at y/s NN = 62.4 GeV scaled by the number of 
valence quarks in the hadron (n q ) and plotted versus pr/riq (a) and (mi — mo)/n q (b). In each case a polynomial curve is fit 
to all particles except pions. The ratio of vifriq to the fit function is shown in the bottom panels (c) and (d). 
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FIG. 6: (color online). V2/n q scaled by the mean eccentricity of the initial overlap region versus (mr — mo)/n, for 0%— 10%, 
10%-40%, and 40%-80% most central Au+Au collisions at v /s^ A 7 = 62.4 GeV. 



plane analysis method has been used. In the bottom 
panels of Fig.[7]we show the ratio v\jv\ for charged pions, 
neutral kaons, and hyperons. The uncertainty in v^jv\ 
from possible non-flow leads to asymmetric errors. The 
ratio V4/V2 is well above 0.5 even when errors are taken 
into account. 

In simple coalescence models (3jj, the ratio v^jv\ for 
hadrons is related to v^jv\ for quarks: 

r , 9-1 Meson , , . r , on Quark , - 

h/« 2 2 ] 2pT w 1/4 + (1/2) [v A /vl] * (2) 



[V-l]™ M W + (V3) (3) 
where here px is the quark pj>. The V^jv\ for mesons can 



also be related to v±jv\ for baryons: 

r , 91 Baryon /rt /n ln * r , 9-1 Meson , 

h/v%] 3p J ^l/6+(2/3)[v 4 /vj] 2pT (4) 

Within this simple model, the large v^jv\ ratios pre- 
sented here indicate a large quark V4. At intermediate 
Pt, where quark-scaling is thought to be valid, we use the 
equations above to fit v±jv\ simultaneously for mesons 
and baryons, with v 4 jv\ for quarks as a free parame- 
ter. The fit range is pr > 1.2 GeV/c for mesons and 
Pt > 1.8 GeV/c for baryons. A good \ 2 P er degree-of- 
freedom (4.4/13) is found with quark u 4 /t>| = 1.93±0.29. 
The grey bars in the bottom panels of Fig.[7]show the cor- 
responding V4/V2 values for mesons and baryons. (v^/v^) 
values for pT/n q > 0.6 GeV/c from data and the fit 
are listed in Table [TTJ Since pion V2 is known to devi- 



10 



■* 0.02 



+ 71* (dE/dx) 
LK° S 

OK* (dE/dx) 




Transverse Momentum p T (GeV/c) 



FIG. 7: (color online). Top panels: minimum bias 114 for pi- 
ons, charged kaons, Kg, anti-protons and A + A at %/s NN — 
62.4 GeV. In the left panel the solid (dashed) line shows the 
value for v\ for pions (kaons). In the right panel the dashed 
line is v\ for A + A. Bottom panels: 114 scaled by v\ (points 
where 114 and V2 fluctuate around zero are not plotted). Grey 
bands correspond to the fit results described in the text and 
Table ITT] The systematic errors on the V4/V2 ratio from non- 
flow are included in the error bars leading to asymmetric er- 
rors. 



ate from the simple scaling laws, we also performed the 
fit excluding the pion data points (fit II) which yielded 
a V4/V2 — 2.18 ± 0.40 and x 2 per degree-of-freedom of 
2.3/9. The small x 2 values for both fits indicate that 
our data are consistent with quark-number scaling where 
quark V4/V2 is approximately 2. 
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(color online). Top panel: V2 for pions and protons 
■ = 62.4 and 17.3 GeV. The 62.4 GeV data are from 



TOF and dE/dx measurements combined. Middle and bot- 
tom panel: ratios of V2 for ty + + n~ , Kg, p + p, A + A and 
at different center-of-mass energies scaled by the values at 
62.4 GeV. The grey and yellow bands represent systematic 
uncertainties in the V2 ratios arising from non-flow effects. 
The grey bands (above unity) are the uncertainties for the 
200 GeV/62.4 GeV data and the yellow bands (below unity) 
are for the 17.3 GeV/62.4 GeV data. 



TABLE II: The ratio w 4 /«l for p T /n q > 0.6 GeV/c from a 
combined fit and from data. Pion data points are used for 
fit I and excluded for fit II. The x 2 P er degree-of-freedom is 
also shown on the bottom row. 
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E. Collision Energy Dependence 

In Fig. [8] (top panel) we plot pion and proton V2 
from = 62.4 Au+Au and 17.3 GeV Pb+Pb col- 

lisions 14]. In the bottom panels we show pion, K 8, 
proton, and A + A data from 17.3 and/or 200 GeV scaled 
by 62.4 GeV data. The 200 to 62.4 GeV ratios are taken 
using v 2 data measured within the 0%-80% centrality 
interval. The TOF v 2 measurements presented in this 



article allow us to show the 17.3 GeV to 62.4 GeV v 2 
ratio to higher px than the 200 GeV data extends. In or- 
der to approximately match the centrality interval used 
for the 17.3 GeV data, the 17.3 to 62.4 GeV ratios are 
taken using respectively 0%-43.5% and 0%-40% central- 
ity intervals. The STAR data at 62.4 and 200 GeV are 
measured within the pseudo-rapidity interval \rj\ < 1.0 
and the 17.3 GeV data are from the rapidity interval 
< y < 0.7. These intervals represent similar y/ybeam 
values. The same method is used to analyze the 200 and 
62.4 GeV data. 

Systematic errors from weak-decay feed-down and 
tracking errors will mostly cancel when taking the ra- 
tio of v 2 at 200 and 62.4 GeV. Possible non-flow errors 
are larger at 200 GeV than at 62.4 GeV. In the lower 
panels of Fig. [51 the shaded bands around unity show 
the uncertainty in the energy dependence of the v 2 ratio 
arising from possible changes in the magnitude of non- 
flow effects at different energies. The portion of the band 
above unity applies to the ratio of 200 and 62.4 GeV data 
while the portion below unity only applies to the ratio of 
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the 17.3 and 62.4 GeV data. 



CONCLUSIONS 



The v 2 data for pions and kaons at 62.4 GeV tends to 
be about 5% smaller than the 200 GeV data (although 
at pt > 1 GeV/c the difference is within systematic un- 
certainties). The anti-proton data at 62.4 and 200 GeV 
are consistent within errors. The data exclude a proton 
v 2 variation between 62.4 and 200 GeV greater than ap- 
proximately 15%. The A + A data show a potentially 
interesting p^-dependence: for px < 1.5 GeV/c the 200 
GeV A + A v 2 is systematically smaller than the 62.4 
GeV data while for p T > 1.5 GeV/c the 200 GeV A + A 
v 2 data are consistent with or larger than the 62.4 GeV 
data. Such a dependence can arise if the system in 200 
GeV collisions develops a larger expansion velocity. 

Appreciable differences are seen between the 17.3 GeV 
and 62.4 GeV data. At p T > 0.5 GeV/c, for both pi- 
ons and protons, the v 2 values measured at 62.4 GeV are 
approximately 10%-25% larger than those measured at 



17.3 GeV [1J, |41(. Although the magnitude of v 2 is dif- 



ferent at the lower energy, the systematics of the particle- 
type dependencies are similar. In particular, pion v 2 and 
proton v 2 cross over each other (or attain similar values) 
at p T near 1.7 GeV/c for = 17.3, 62.4 and 200 GeV 

data. Due to the limited kinematic range covered by the 

17.3 GeV data, a quark-number dependence of v 2 at in- 
termediate pt can neither be confirmed nor excluded. 

The increase in the magnitude of v 2 from 17.3 GeV to 

62.4 GeV and the similarity of 62.4 GeV v 2 to 200 GeV 
v 2 has been taken as a possible indication for the onset 
of a limiting behavior [42|]. In a collisional picture, a sat- 
uration of v 2 could indicate that for ^/s NN at and above 
62.4 GeV the number of collisions the system constituents 
experience in a given time scale can be considered large 
and that hydrodynamic equations can therefore be ap- 
plied. Hydrodynamic model calculations of v 2 depend 
on the model initialization and the poorly understood 
freeze-out assumptions [13, [Tlj] . As such, rather than 
comparing the predicted and measured values at one en- 
ergy, the most convincing way to demonstrate that a hy- 
drodynamic limit has been reached may be to observe 
the onset of limiting behavior with y/s NN . For this rea- 
son, v 2 measurements at a variety of center-of-mass en- 
ergies are of interest. Contrary to the large differences 
reported in Ref. [i^], we find that when the 17.3 and 62.4 
GeV v 2 (pt) data are compared within similar \y\/ybeam 
ranges and when possible non-flow systematic uncertain- 
ties are accounted for (the yellow bands in the bottom 
panel of Fig. [8]), the differences between v 2 (pt) within 
the data sets may be as small as 10%-15%. As such, 
a large fraction of the deviation between the SPS data 
and hydrodynamic models arises due to the wide rapid- 
ity range covered by those measurements (v 2 approaches 
zero as beam rapidity is approached |43j), increased (pt) 
values at RHIC and the larger v 2 values predicted for the 
lower colliding energy by hydrodynamic models. 



We presented measurements of v 2 for pions, kaons, pro- 
tons, A, A, S + S, and f2 + fi from Au+Au collisions with 
ys^" = 62.4 GeV. We compared these measurements 
to similar measurements at ^s NN — 17.3 and 200 GeV. 
The 62.4 GeV pion, kaon, proton, and hyperon v 2 data 
are, within a few percent, consistent with the equivalent 
data at 200 GeV. Within similar y/ybeam intervals and 
and after we account for systematic uncertainties, we find 
that for a given identified particle species the difference 
between 17.3 and 62.4 GeV v 2 data may be as small as 
10%-15%. We find that A v 2 is larger than A v 2 at 62.4 
and 200 GeV and that the difference is larger at 62.4 GeV 
where the anti-baryon to baryon yield ratio is smaller. At 
both energies our measurements are consistent with net 
A v 2 being approximately 10%-15% larger than A and 
pair-produced A v 2 . 

Our v 2 measurements at 62.4 GeV are consistent 
with the quark-number scaling of v 2 first observed from 
Au+Au collisions at ^/s^; = 200 GeV. The 17.3 GeV 
data do not extend to high enough pt to test quark- 
number scaling. We note, however, that the pt where 
the v 2 values for mesons and baryons cross over each 
other (or, in the case of 17.3 GeV data, become similar) 
is approximately the same at all three center-of-mass en- 
ergies. This indicates that identified particle v 2 at 17.3 
GeV may also be consistent with quark-number scaling. 

We also reported measurements of the higher harmonic 
term, ^4, for pions, kaons, protons, and A+A. These mea- 
surements are also consistent with quark-number scal- 
ing laws arising from coalescence or recombination mod- 
els [39j . This quark- number dependence may indicate 
that in ultra-relativistic heavy-ion collisions collective 
motion is established amongst quarks and gluons be- 
fore hadrons are formed. This view is supported by 
the large v 2 values measured for multi-strange baryons 
at v /s^7 = 62.4 and 200 GeV [21]. Collisions involv- 
ing lighter nuclei and larger, deformed nuclei (U+U) will 
provide another opportunity to study mass and quark 
number systematics for v 2 . The possible approach to 
limiting values for v 2 (where the pt and mass depen- 
dence at pt < 1 GeV/c are consistent with hydrody- 
namic models) along with the evidence presented here 
that the relevant degrees of freedom in the early sys- 
tem may be sub-hadronic (e.g. constituent quarks) sug- 
gests that a strongly coupled matter with sub-hadronic 
degrees of freedom may be created in heavy-ion collisions 
at RHIC. 
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